When diffraction rules the stereodynamics of rotationally inelastic collisions 
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Following upon our recent work on vector correlations in the Ar-NO collisions [PCCP 12, 1038 
(2010)], we compare model results with close-coupling calculations for a range of channels and colli- 
sion energies for the He-NO system. The striking agreement between the model and exact polariza- 
tion moments indicates that the stereodynamics of rotationally inelastic atom-molecule collisions at 
thermal energies is governed by diffraction of matter waves from a two-dimensional repulsive core 
of the atom-molecule potential. Furthermore, the model polarization moments characterizing the 
He-NO, He-C>2, He-OH, and He-CaH stereodynamics are found to coalesce into a single, distinctive 
pattern, which can serve as a "fingerprint" to identify diffraction-driven stereodynamics in future 
work. 
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The pioneering work of Herschbach and coworkers [T] 
on vector correlations in the domain of molecular colli- 
sions spurred an effort to extract directional information 
hidden in molecular dynamics experiments and compu- 
tations, and thereby to enrich our knowledge of how a 
given collision proceeds [5] . However, even when charac- 
terized to the full by vector correlations, the why of col- 
lision dynamics can only be answered as well as the the- 
oretical method applied to treat the collision allows [3]. 
Therefore, we implemented an analytic model of colli- 
sion dynamics, capable of answering the why for a class 
of collisions in detail, and used it to develop an ana- 
lytic model of vector correlations in such collisions [I]. 
The collision model is based on the Fraunhofer scatter- 
ing of matter waves [5], recently extended to treat colli- 
sions in fields [5] as well as the stereodynamics of rota- 
tionally inelastic atom-molecule collisions at thermal and 
hyperthermal energies [1]. The Fraunhofer model [1HS] 
relies on the sudden approximation, which treats the ro- 
tational motion as frozen during the collision and thereby 
allows to replace the inelastic scattering amplitude with 
the elastic one. The elastic scattering amplitude itself 
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is approximated by the amplitude for Fraunhofer diffrac- 
tion of matter waves from a sharp-edged, impenetrable 
obstacle acting in place of the molecular scatterer and 
captures forward scattering. At collision energies of hun- 
dreds of cm , consistent with the sudden approxima- 
tion, the shape of the scatterer is approximated by the 
repulsive core of the atom-molecule potential, with the 
attractive part disregarded. The Fraunhofer model ren- 
ders fully state- and energy-resolved scattering ampli- 
tudes and all the quantities that unfold from them in 
analytic form. 

The vector correlations obtained from the Fraun- 
hofer model were found to closely reproduce those ex- 
tracted from close-coupling calculations and experiment 
of Wade et al. [7] for the much examined Ar-NO (A" 2 11) 
system [9]. This agreement allowed interpreting the 
system's collision stereodynamics in terms of the Fraun- 
hofer model. The Fraunhofer model of vector correlations 
henceforth revealed that the stereodynamics of the Ar- 
NO rotationally inelastic collisions is dominated by the 
diffractive part of the scattering amplitude which is gov- 
erned by a single Legendre moment characterizing the 
anisotropy of the hard-core part of the system's poten- 
tial energy surface. Given the "geometric" origin of this 
behavior - ordained by the angular momentum disposal 
- we wondered about its generality. In this work, we 
compare model results with close-coupling calculations 
for different channels and collision energies of the He- 
NO system. The striking agreement between the Fraun- 
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hofer model and the close coupling calculations found 
herein and in our previous work on the Ar-NO collision 
stereodynamics [3] attests to the dominant role of diffrac- 
tion in rotationally inelastic collisions. Furthermore, the 
model alignment moments were found to exhibit remark- 
able similarities for different collision partners, such as 
He-NO, He-0 2 , He-OH, and He-CaH, which allowed us 
to identify the form factors ( "fingerprints" ) of diffraction- 
driven stereodynamics. 

We note that the discrepancy between the model and 
the exact calculations of the differential cross sections 
arises from the non-diffractive contributions to scattering 
rather than from the neglected diffraction by the long- 
range potential @]. 

The stereodynamics of an atom-diatom collision is 
usually described by a set of four vectors: the initial and 
final relative velocities, k and k', and the initial and final 
rotational angular momenta of the diatomic molecule, j 
and j'. We use the initial and final relative velocities k 
and k' to define the XZ plane of the space-fixed coor- 
dinate system, with the initial relative velocity k point- 
ing along the Z axis. In keeping with the convention of 
de Miranda et al. [10) . we characterize the spatial dis- 
tribution of the angular momenta relative to the XYZ 
frame by the real polarization moments a,q±, which are 
related to the multipole moments in the expansion of 
the density operator in terms of the irreducible tensor 
operators [llj . Since, within the Fraunhofer model, the 
scatterer is two-dimensional, the model can only account 
for even-fc (alignment) polarization moments with even 
q [H |6] . The "extra symmetry" of the model causes odd- 
k (orientation) moments as well as odd-q alignment mo- 
ments to vanish. 

First, we focus on vector correlations in the He-NO 
system, for different scattering channels and collision en- 
ergies. We compare analytic model results with fully 
quantum close-coupling calculations performed using the 
HIBRIDON suite of computer codes [12] on the PES of 
Klos et al. [13) . In order to characterize the k — k' — j' 
three-vector correlation, we make use of the alignment 
moments Oq(j') and of the diatomic's final rota- 

tional angular momentum j' with respect to the XYZ 
frame. The physical meaning and range of the polariza- 
tion moments is given in Table [I] 

Figure [l] shows the ag(j') and a| + (j') moments for dif- 
ferent channels of the He-NO (j = Q = 1/2 -> /, Q' = 
1/2) system at a collision energy of 520 cm -1 . All the 
alignment moments we present below were obtained for 
unresolved initial and final lambda-doublet states. The 
compelling agreement of the Fraunhofer model with the 
exact calculations for the He-NO as well as for Ar-NO [J] 
systems attests to the predominant role of diffraction in 
shaping the stereodynamics of rotationally inelastic col- 
lisions at thermal and hyperthermal energies. Further- 
more, the agreement with the model shows that what 
matters the most - as far as the PES is concerned - is 
the 2D contour of its repulsive core. Save for the mo- 
ments' oscillatory structure, which differs for different 



TABLE I: The physical meaning and range of the ao(j') and 
a 2+(j') alignment polarization moments. The Z axis points 
along the initial relative velocity k. The final relative veloc- 
ity k' lies in the X > half of the XZ plane. The ag(j") 
moment accounts for alignment of j' with respect to k. Pos- 
itive (negative) values of Oo(j') correspond to j' || k (j' _L k, 
in which case j' lies in the XY plane). The positive (nega- 
tive) values of the ffll+0') moment correspond to alignment 
of j' along the X-axis (F-axis). The a2_(j') moment vanishes 
identically. The indicated ranges of the moments correspond 
to the high-j limit. 



Moment 


ao(j') 


«2 + (j') 


Meaning 
Range 


j' along Z 

j'±Z -> -1/2 
i'\\Z -> 1 


j' along Y or X 

j'\\Y -> -y/3/2 
j'\\X -> VS/2 



final states (as these arise due to different sets of the 
PES's Legendre moments), the diffraction manifests it- 
self in the same way in all the scattering channels, i.e., 
it leaves behind the same fingerprints, shown in Fig. [I] 
by the thick semi-transparent lines: whereas the Oq mo- 
ments are negative for a purely forward scattering and 
increase at larger scattering angles where they approach 
a constant positive value, the a\ + moments drop from a 
zero for forward scattering and tend to a constant nega- 
tive value at larger scattering angles. This means that, 
due to geometric reasons, j' aligns perpendicular to k for 
a purely forward scattering. At larger scattering angles, 
j' aligns along the Z and Y axes, indicating that the 
molecular axis aligns preferentially in the X direction. 

If the oscillations of the alignment moments are due to 
diffraction, they should scale with the de Broglie wave- 
length of the collision system and the size of the molec- 
ular scatterer, in analogy with the wavelength of light 
and the obstacle size in optics. Indeed, Figure [2] reveals 
such a scaling of the ag(j') and moments for the 

He-NO {j = fl = 1/2 f = 9/2,0' = 1/2) channel 
and a range of collision energies. The scaling was im- 
plemented by the transformation $ —¥ i9i?o/A, where $ is 
the scattering angle, Rq is the molecular size and A is the 
de Broglie wavelength. Whereas the alignment moments 
furnished by the Fraunhofer model coincide exactly upon 
scaling in ■§ (black solid line) , such a scaling brings the ex- 
act moments (filled circles, squares, and triangles) quite 
close to one another, but does not result in their exact 
matching. The remaining differences among the scaled 
exact moments corresponding to different collision ener- 
gies attest to non-diffractive contributions to scattering, 
the influence of the attractive branch of the PES, and a 
breakdown of the sudden approximation. For instance, 
the scaled exact moments corresponding to 520 and 1000 
cnr 1 come close to one another, while the oscillations for 
150 cm -1 are quite off, especially for the czqQ') moment. 
We ascribe this discrepancy to the potential well of 25 
cm -1 [13j that brings about non-diffractive contributions 
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FIG. 1: Polarization moments pertaining to the He-NO (j = fi = 1/2 — > j' , fi' = 1/2) system at 520 cm -1 . Model results are 
shown by lines, exact computations - by symbols. Form factors of the alignment moments are shown by thick semitransparent 
lines. 
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FIG. 2: Polarization moments versus scaled scattering angle 
pertaining to the He-NO (j = Q, = 1/2 -> j' = 9/2, fi' = 1/2) 
system at different collision energies. The Fraunhofer mo- 
ments coincide exactly. Form factors of the alignment mo- 
ments are shown by thick semitransparent lines. See text. 



to the stereodynamics while, at the same time, diminish- 
ing the role of the PES's "repulsive core." However, the 
form factor of the alignment moments (thick semitrans- 
parent line) is quite similar to the ones shown in Fig- 
ure [T] indicating that diffraction leaves behind the same 
"fingerprint" for different collision energies and scatter- 
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FIG. 3: Polarization moments versus scaled scattering angle 
pertaining to the He-NO (j = f2 = 1/2 -> f = 9/2, SI' = 
1/2), He-0 2 (j = 0,N = l,-> f = A,N' = 5), He-OH (j = 
Q = 1/2, -> j' = 9/2, fi' = 1/2), and He-CaH (j = 1/2, N = 
0, -> j' = 11/2, N' = 6) systems at 520 cm" 1 . The moments 
were obtained using the Fraunhofer model. Form factors of 
the alignment moments are shown by thick semitransparent 
lines. See text. 



ing channels. Interestingly, at low collision energies (10 
cm" 1 ), the oscillations of the alignment moments fall fur- 
ther out of phase, but the form factors still remain in 
place. 
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In order to see how the diffraction patterns change 
from one scattering system to another, we also examined 
the He-C>2, He-OH, and He-CaH systems at a collision 
energy of 520 cm -1 and similar channels, using potential 
energy surfaces of refs. . 

Figure [3] shows the dependence of the polarization mo- 
ments on the scaled scattering angle obtained from the 
Fraunhofer model, pertaining to the He-NO (j = O = 
1/2 = 9/2, Q' = 1/2), He-0 2 (j = 0, N = 1, ->• f = 
4,7V' = 5), He-OH ( j = ft = 1/2, -s- f = 9/2,17' = 1/2), 
and Hc-CaH (j = 1/2, N = 0,-> j' = 11/2, AT' = 6) 
systems. One can see that the shape of the oscillations 
differs for different collision partners. Going from the 
most symmetric scatterer, O2, through NO and OH to 
the most asymmetric one, CaH, one can see that the os- 
cillations become increasingly asymmetric too, and their 
amplitudes decrease. This suggests relating the asymme- 
try and amplitude of the oscillations to the asymmetry 
of the repulsive core of the PES. However, the form fac- 
tors (shown by the semi-transparent curves) that capture 
the alignment moments are very similar to one another, 
indicating that the "fingerprints" of diffraction barely de- 
pend on the collision system. We note that the same can 
be said about the higher-rank alignment moments, which 
we evaluated as well [T7] . 

In summary, in order to gain insight into the stere- 



odynamics of rotationally inelastic atom-molecule colli- 
sions, we compared polarization moments obtained from 
an analytic model with those extracted from exact close- 
coupling calculations for the He-NO collisions. The 
model alignment moments were found to come as close 
to exact results as for the previously examined Ar-NO 
system , which reveals that the collision stereodynam- 
ics in question is governed by diffraction of matter waves 
from a 2D contour of the repulsive core of the potential 
(flat sharp-edged obstacle). The oscillatory patterns of 
the alignment moments due to diffraction scale with the 
de Broglie wavelength and the molecular size. Therefore, 
deviations from such patterns single out other contribu- 
tions to the scattering which are mainly non-diffractive. 
Furthermore, diffraction leaves behind the same finger- 
prints for different channels and collision energies for 
a range of systems, including He-NO, He-02, He-OH, 
and He-CaH. These fingerprints can be used to identify 
diffraction-driven stereodynamics in future experiments 
and exact computations. 
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